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ABSTRACT 

It is shown that collisionless shock waves can be driven in unmagnetized electron-positron plasmas 
by performing a two-dimensional particle-in-cell simulation. At the shock transition region, strong 
magnetic fields are generated by a Weibel-like instability. The generated magnetic fields are strong 
enough to deflect the incoming particles from upstream of the shock at a large angle and provide an 
effective dissipation mechanism for the shock. The structure of the collisionless shock propagates at 
an almost constant speed. There is no linear wave corresponding to the shock wave and therefore this 
can be regarded as a kind of "instability-driven" shock wave. The generated magnetic fields rapidly 
decay in the downstream region. It is also observed that a fraction of the thcrmalizcd particles in 
the downstream region return upstream through the shock transition region. These particles interact 
with the upstream incoming particles and cause the generation of charge-separated current filaments 
in the upstream of the shock as well as the electrostatic beam instability. As a result, electric and 
magnetic fields are generated even upstream of the shock transition region. No efficient acceleration 
processes of particles were observed in our simulation. 

Subject headings: shock waves — plasmas — instabilities — magnetic fields — acceleration of particles 



1. INTRODUCTION 

Collisionless shocks are driven in various situations 
in astrophysical plasmas. They often accelerate par- 
ticles and generate nonthermal high-energy particles. 
Such shocks in supernova remnants are known as 
the accelerators of n onthermal electrons (e.g., SN1006; 
iKovama et al. Ifl995[) . and are also believed to be the ac- 
celerators of cosmic rays with energies up to the knee 
energy (~ 10 15 5 eV). Collisionless shocks also exist in 
electron-positron plasmas, for example, that in the Crab 
nebula, and they are considered to generate nonthermal, 
high-energy electrons and positrons as well. 

Recently, it has been suggested from observations 
that the magnetic fields are amplified or generated 
around collisionless sho c ks in several supernova remnants 
(IVink fc Laming 1120031 iBamba et al. 1 120031 : 1 Yolk et al. I 
|2005|) . These magnetic fields may be generated b y the 
high- energy particles accelerated at the shocks fe.g.. lBelfl 
I2004D . On the other hand, such a generation mechanism 
can be related to the microphysics of the collisionless 
shocks themselves. For the generation of magnetic fields 
in the relativistic s hocks associated with the af terglows of 
gamma-ray bursts, Mcdvcd ev fc Loebl (Il999h suggest ed 
that the Weibel instability (jWeibel lll959HFried Ill959f) is 
driven at the shock and generates strong magnetic fields. 
This mechanism can work in shocks in electron-positron 
plasmas as well. 

For electron-positron plasmas, the Weibel instabil- 
ity has been investigated in detail by using numerical 
simulations. Several particle-in-cell simulations showed 
that in counterstreaming electron-positron plasmas, the 
Weibel instability develops an d generates a strong 
sub-equipartition magnetic field (jKazimura et al. 1 1 19981 : 
iSilva et al. 1 l2003h . Furthermore, some authors found 
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that in their simulations with longer calculation times, 
a shock-like structure that associates with a strong mag- 
netic field is forme d even in unma g netized electron- 
posit ron plasmas (jGruzinov ] 120011 : lHaruki fc Sakail 
2003). Thus, it should be clarified whether the shock-like 
structure can be regarded as a shock from a macroscopic 
point of view, for example, whether the structure dis- 
sipates the upstream bulk kinetic energy at the "shock" 
transition region or propagates at a steady speed into the 
upstream plasma, etc., as well as other detailed features. 
If they are considered as shocks, it is also important to 
elucidate the role of the strong magnetic fields generated 
by the Weibel-like instability in the dissipation of the 
shock. 

In this study, we further investigate the dynamics 
of collisionless shocks in unmagnetized electron-positron 
plasmas in detail by performing a high-resolution two- 
dimensional particle-in-cell simulation and confirm that 
a kind of collisionless shocks indeed forms mainly due to 
the magnetic fields generated by the Weibel-like instabil- 
ity. However, we also find that no efficient particle accel- 
eration occurs in the shock and the generated magnetic 
fields rapidly decay within the shock transition region, 
while in the real world, particles are accelerated at the 
shocks and large-scale, sub-equipartition magnetic fields 
are generated around the shocks. The reasons why these 
processes are not observed in our simulation may be re- 
lated to the background magnetic field, the spatial and 
temporal scales of the simulation, or the composition of 
the plasma. 

2. SIMULATION 

In order to investigate the collisionless shocks in 
electron-positron plasmas without background magnetic 
fields, we performed numerical simulations. The simu- 
lation code is a relativistic, electromagnetic, particle-in- 
cell code with two spatial and three velocity dimensions 
(2D3V) developed based on a standard method described 
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bv iBirdsall fc Langdonl (|1991( ). Thus, the basic equa- 
tions of the simulation are the Maxwell's equations (in 
Gaussian units): 

10E „ 4tt, ldB 

-— = VxB J, - =-VxE, 1 

c at c c at 

with 

V • E = 47rp, V • B = 0, (2) 

as constraints, where c is the speed of light; E and B, 
the electric and magnetic fields, respectively; J, the cur- 
rent density; and p, the charge density, together with the 
equation of motion of particles: 

at mc \ 1 J 

where q and m are the charge and mass of the particle, 
respectively; 7, the Lorentz factor of the particle; and 
u = 7v/c, the 4- velocity of the particle (where v is the 
ordinary 3- velocity) . Particle velocities are often denoted 
in terms of the 4-velocity throughout this paper. 
In the following, we consider t = ui~ e0 as the unit 

of time and the electron skin depth Iq = cu;~ as the 

unit of length, where a; po ,o = (47rn e oe 2 /m e ) 1 / 2 is the 
electron plasma frequency defined for the mean electron 
number density n e o; m e is the electron mass, and — e is 
the electron charge. The units of electric and magnetic 
fields are both taken as E* = £?* = c(47rn e om e ) 1//2 ; they 
are defined so that their corresponding energy densities 
are both equivalent to the half of the mean electron rest 
mass energy density. 

The simulations were performed on a grid of N x x N y = 
4096 x 512 with 3.2 x 10 8 particles for each species. The 
physical size of the simulation box is L x x L y = 480 x 60 
in the unit of Iq. Thus, the size of a cell is (Ax, Ay) = 
(0.12,0.12) in the same unit. The time step is taken as 
At = 0.025. 

Since we consider collisionless shocks in unmagnetized 
plasmas, the electromagnetic fields were initially set to 
zero over the entire simulation box. The boundary con- 
dition of the electromagnetic field is periodic in each di- 
rection. 

In this simulation, a collisionless shock wave is driven 
according to the well-known "injection method." There 
are two walls at x — 30 and 450. They reflect parti- 
cles specularly, but let the electromagnetic waves pass 
through freely; the electromagnetic waves propagating 
outside the walls (0 < x < 30 and 450 < x < 480) 
are dissipated by means of Ohmic dissipation so that 
they cannot reach the wall at the other side. Initially, 
both the electrons and positrons are loaded uniformly in 
the region between the two walls with a bulk velocity of 
u x = 2.0 (thus, the corresponding bulk Lorentz factor is 
r = 2.24). Their thermal velocities were given so that 
the distribution of each component of the 4-velocity Ui 
(i = x, y, z) obeys the Gaussian distribution with a stan- 
dard deviation of crth = 0.1 in the plasma rest frame, or 
the "fluid" frame. (In effect, this generates the Maxwell 
distribution.) At the initial stage of the simulation, par- 
ticles that were located near the right wall (at x = 450) 
were reflected by the wall and then interacted with the 
incoming particles, i.e., the upstream plasma. This in- 
teraction causes some instabilities and then a collision- 
less shock is formed. It should be noted that the frame 
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Fig. 1. — Electron number density at t = 233. The value is 
normalized to the far upstream number density measured in the 
downstream rest frame (i.e., the simulation frame). The left-hand 
side and the right-hand side of this figure correspond to upstream 
and downstream of the shock, respectively. The transition region 
extends from x = 350 to 400. 




Fig. 2. — Time development of the electron number density. 
The horizontal axis is the x coordinate and the vertical axis is the 
time. The plotted electron number density is averaged over the 
y direction. [See the electronic edition of the Journal for a color 
version of this figure.] 

of the simulation corresponds to the downstream rest 
frame. Thus, we observed the propagation of the shock 
from right to the left in the downstream rest frame. 

3. RESULTS 

The following figures show the quantities at the end of 
the simulation (t = 233) otherwise stated. 

3.1. Basic Features of Collisionless Shock 

Figure [1] shows the electron number density normalized 
to the initial number density, or the far upstream num- 
ber density, measured in the downstream rest frame, n e Q. 
In this figure, the left-hand side and the right-hand side 
correspond to upstream and downstream of the shock, re- 
spectively. We observe that the number density rapidly 
increases while crossing the transition region, which ex- 
tends roughly from x — 350 to 400, and then it becomes 
homogeneous by moving further downstream. The tran- 
sition region does not have a one-dimensional structure 
but a two-dimensional filamentary structure with density 
fluctuations along the y direction. The typical size of the 
filaments is in the order of several electron skin depths. 

Figure [2] shows the time development of the electron 
number density. The horizontal and vertical axes repre- 
sent the x coordinate and time, respectively. The plotted 
number density is averaged over the y direction. We ob- 
serve that the transition region, which is visible as the 
jump in the number density, propagates upstream, i.e., 
to the left, with time at an almost constant speed. From 
this figure, the propagation speed measured in the down- 
stream rest frame is estimated to be V^h.d ~ —0.39c. 

The thin structure propagating upstream at almost the 
speed of light, which is visible in the lower-left portion 
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Fig. 3. — The y-averaged profiles of quantities at t = 233: (a) 
the electron number density n e /n eu (solid curve for t = 233 and 
dashed curve for t = 200), (b) the mean velocity in the a;-direction 
V x /c, and (c) the magnetic field \B\/B* (solid curve) and electric 
field \E\/E* (dashed curve). 

of this figure, is due to the particles reflected at the right 
wall at a very early stage of the simulation. This is of 
course the consequence of the initial condition, but such 
a situation may be realized in real cases. Anyway, this 
structure fades away with time, and it would not affect 
the structure around the transition region later. 

Figure [3] shows the y-averaged profiles of three quanti- 
ties around the shock transition region: (a) the electron 
number density, n e /n e o; (b) the mean velocity of elec- 
trons in the x-direction, V x /c] and (c) the magnitude 
of electric and magnetic fields, \E\/E* and \B\/B*. In 
Fig. [3t, the solid curve and the dashed curve represent 
the number densities at t = 233 and t = 200, respectively. 
Comparing the two curves, we see that the transition re- 
gion propagates upstream keeping its averaged structure 
almost unchanged. The compression ratio measured in 
the downstream rest frame is approximately 3.3. The 
mean velocity of electrons in the x direction, V Xl rapidly 
approaches zero within the transition region (Fig.[3jD). In 
Fig. [3t, it is notable that a strong magnetic field exists 
within the transition region (350 < x < 400). We also 
see that behind the shock transition region, the strength 
of the magnetic field decreases and strong electric and 
magnetic fields exist even in the upstream region. These 
points will be mentioned later. 

FigureSt shows the z component of the magnetic field, 
B z . It is clear that a strong magnetic field exists within 
the transition region (350 < x < 400), as indicated in 
Fig. St- Figure HJd shows the current density in the x 
direction, J x . We observed that there are considerable 
number of current filaments that carry currents in the 
x direction within the transition region. The filamen- 
tary structure observed in the number density (Fig. [T]) 
is in fact due to the existence of these current filaments. 




Fig. 4. — (a) The z component of the magnetic fields B z and 
(b) the x component of the current density J x at t = 233. 



They generate a magnetic field fluctuating in the y di- 
rection, as shown in Fig. 2t. This situation in which 
strong magnetic fields fluctuating perpendicular to the 
direction of streaming are generated is similar to those 
in the simulations o f the Weibel instability in counter- 
strea ming plasmas (|Kazimura et al. I Il998t ISilva et aJ~l 
2003). Therefore, it is reasonable to consider that in 
the transition region in which the upstream plasma is 
mixed up with the downstream plasma, a large velocity 
dispersion in the x direction and hence a large anisotropy 
in the particle velocity distribution are induced; further, 
the anisotropy drives the Weibel-type instability, which 
generates the magnetic field. The typical magnitude of 
the magnetic field in this region is \B\ ~ 0.7-B* and the 
corresponding Larmor radius for the incoming upstream 
particles with u — 2.0 is r g ~ 2.9, which is comparable to 
the size of the magnetic field structure, i.e., the current 
filaments, in the transition region. Therefore, the in- 
coming upstream particles should be isotropised in this 
region because they are deflected at a large angle by the 
strong magnetic field at this location. This isotropisa- 
tion provides an effective dissipation mechanism for the 
upstream bulk kinetic energy. 

Figure [5] shows the y-averaged profile of the magnetic 
energy density normalized by the upstream bulk kinetic 
energy density. We see that the magnetic energy density 
reaches about 10% of the upstream bulk kinetic energy 
density at the peak and decays downstream of the shock. 
The saturation and decay mechanism of the magnetic 
field would b e similar to those of the Weibel instability 
(|Kato I [20051 ). The current and the magnetic field gen- 
erated by the instability in each filament saturates when 
the magnetic field becomes strong enough to deflect the 
current-carrying particles in the filament. This predicts 
the sub-equipartition magnetic field at saturation and it 
is consistent with the simulation result. After the satura- 
tion, the current filaments will coalesce with each other 
to form larger filaments with the decreasing magnetic 
field strength. We observe the evolution of the structure 
to a larger scale together with the decay of the magnetic 
field strength downstream of the shock transition region 
in Figs. Skint, andU 

Figure [6] shows (a) the fluid velocity of electrons in 
the x-direction, V x /c, and (b) the mean kinetic energy 
of the electrons measured in the local fluid rest frame of 
the electrons expressed in terms of electron rest mass, 
(7 — f) fluid- It is evident that the incoming plasma from 
upstream is rapidly decelerated within the transition re- 
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Fig. 5. — The y-averaged profile of the magnetic energy den- 
sity normalized by the upstream bulk kinetic energy density at 
t = 233, eg /rn e om e c 2 . The magnetic field strength reaches a sub- 
equipartition level at the peak and rapidly decays in the down- 
stream region. 
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Fig. 6. — (a) The x component of the electron mean velocity 
V x /c and (b) the mean kinetic energy of the electrons measured in 
the local plasma rest frame of the electrons, (7 — l)fl u id, at t = 233. 
It is evident that the bulk kinetic energy of the upstream plasma 
is converted into the internal kinetic energy within the transition 
region. [See the electronic edition of the Journal for a color version 
of this figure] 

gion, where the strong magnetic field exists (see also 
Fig. [4ji), and its bulk kinetic energy is converted into 
random or "thermal" kinetic energy in the downstream 
rest frame. 

Figure [7] shows the standard deviations of (a) u x of 
electrons, cr x , and (b) u y of electrons, o~ y , both measured 
in the local plasma rest frame of the electrons. We ob- 
serve that the particles are isotropised in the transition 
region due to the strong magnetic field existing there. 
It should be noted that, in this two-dimensional simula- 
tion, particles are isotrop ised only on the x — y p lane and 
not in the z direction (c.f. lHaruki fc Sakai 112003) because 
only the Weibel modes with wave vectors on the x — y 
plane can develop and hence only B z is generated. In 
a three-dimensional simulation, the Weibel modes with 
wave vectors in the z direction should develop equally 
and particles are isotropised in all directions. 

In a macroscopic view, the features shown in the above 
figures would be observed as those of shock waves; the 
shock waves dissipate the upstream bulk kinetic energy 
in the shock transition region and propagate at an al- 
most constant speed into the upstream plasma. In this 
sense, we can regard this as a collisionless shock. In this 
shock, the magnetic field in the transition region gener- 
ated by the Weibel-like instability plays an essential role 
in the dissipation process of the collisionless shock in un- 
magnetized plasmas. It should be noted that the shock 
has no corresponding fundamental linear wave, which is 
required to define the Mach number, in contrast to, for 
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Fig. 7. — Standard deviation of (a) u x and (b) u y of electrons 
both measured in the local plasma rest frame of the electrons. The 
electrons are isotropised when they pass through the transition 
region. [See the electronic edition of the Journal for a color version 
of this figure.] 

example, the perpendicular shocks in magnetized plas- 
mas, in which the magnetosonic wave is the fundamental 
linear wave of the shock. This point is clear from the fact 
that the wave vector of the magnetic field in the transi- 
tion region is in the y direction and not in the x direction 
or the shock normal in the present case and hence it has 
an essentially a two-dimensional structure. Thus, it can 
be said that this is a class of "instability-driven" shock 
waves. For such shocks, the concept of the Mach number 
would not be applicable, and the propagation speed of 
the shock and the width of the transition region would 
be essentially determined by the nature of the instabil- 
ity that provides the dissipation mechanism at the shock 
transition region. 

3.2. Charge-separated Current Filaments and Backward 
Flowing Particles 

From Fig. [4}d, we observe that there are current fila- 
ments even in the upstream region (x < 350). Figure 
shows the charge density p. Comparing these two fig- 
ures, it is evident that these filaments have a net space 
charge. Since in each filament, J x and p have the same 
sign, these filaments are mainly composed of the incom- 
ing upstream particles with u x > 0. Figure [BJa shows the 
y component of the electric field, E y . We observe that E y 
is generated between these filaments due to their space 
charge. As shown in Fig. Hk, these filaments also gener- 
ate magnetic field B z around them. Thus, rather strong 
electric and magnetic fields exist in the upstream region, 
as shown in Fig. |3j;. Figure[8j: shows the x component of 
the electric field. We observed that strong electric fields 
exist at the upstream edge of the transition region and 
that there is a coherent structure in E x in the upstream 
of the transition region (x < 350). 

The region where the charge-separated current fila- 
ments exist would correspond to "the charge s eparation 
layer" suggested bv iMilosavlievic et al. I (|2006j ). As they 
pointed out, these filaments need a small fraction of the 
particles flowing backward to stabilize themselves. Fig- 
ure [HI shows the x — u x phase-space plot of electrons at 
t = 233. We observe that the particles flowing back- 
ward (u x < 0) are supplied around the shock transition 
region (x ~ 350). These particles would stabilize the 
charge-separated current filaments in the charge separa- 
tion layer. They can also drive the electrostatic beam 
instability. The coherent electrostatic field E x upstream 
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Fig. 8. — (a) Charge density p, (b) the y component of the 
electric field E y and (c) the x component of the electric field E x 
at t = 233. 
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Fig. 9. — The x — u x phase-space plot of electrons at t = 233. 
Upstream of the shock transition region (x < 350), a fraction of 
particles flow backward, which causes two-stream type instabilities. 
In the downstream region (x > 370), electrons are thermalizcd 
(along the x direction). 

of the transition region (x < 350; see Fig. [Sfc) would be 
induced by this instability. These particles should also 
contribute to drive the Weibel-type instability around 
the transition region. 

Figure |TD] shows the time development of the number 
density of electrons again as in Fig. [5J but with contours. 
We observe that there is a "precursor" region in which 
the number density is within a range of 1.1 < n e /n e Q < 
1.9 in front of the transition region. In this region, the 
increase in the number density from that in the upstream 
plasma n e o would be mainly due to the electrons flowing 
backward (see also Fig. The width of the precursor 
region is almost constant (Ax ~ 30) after t ~ 100. 

As already shown in Fig. [4)3, there are current fila- 
ments in the transition region. However, they are charge- 
separated only at the filament edges except those at 
the downstream side that are continuously connected to 
the thermalized downstream region (see Fig. [8^). The 
strength of the magnetic and electric fields, current den- 
sity, and charge density all attain their maximum values 
in this region, and the upstream kinetic energy is dis- 
sipated mainly by the magnetic field and partly by the 
electric field. 

A fraction of the downstream particles flowing back- 
ward with u x < u s h,d ~ —0.42, which is the 4-velocity 
corresponding to V^h.dj can return upstream through the 
filaments that carry the currents in the same direction 
as the particles. On the other hand, if the particles re- 




FlG. 10. — Time development of the electron number density 
as in Fig. [2] but with contours. The width of the precursor region 
(1.1 < n e /n e o < 1-9), where the density of the electrons flowing 
backward is high, is almost constant (Ax ~ 30) after t ~ 100. [See 
the electronic edition of the Journal for a color version of this 
figure.] 
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Fig. 11. — Distribution of backward flowing electrons with 
u x < —2. They move upstream through the current filaments with 
J x > 0. This is evident for x > 350. 

turn through the filament carrying opposite currents, the 
magnetic fields in the filaments repel the particles from 
the filaments. Figure [TT] shows the distribution of elec- 
trons that are flowing upstream with u x < —2. By com- 
paring with J x in Fig. [4}}, we observed that these elec- 
trons indeed flow in the current filaments with J x > 0. 
In a similar manner, the positrons flow backward in the 
current filaments, but with J x < (not shown here). 
However, only a small fraction of the backward flowing 
particles can move to the upstream region (x < 350) be- 
cause large potential barriers exist, caused by the charge 
density at the upstream edge of the filaments (see also 
Figs. EH and[8t). The particles that could pass through 
to the upstream region would contribute to stabilize the 
charge-separated current filaments in the charge separa- 
tion layer, as mentioned earlier. 

3.3. Particle Energy Distribution and Acceleration 

Figure [12] shows the energy histogram of the elec- 
trons accumulated over the area between x — 430 and 
X = 442 in the downstream region (the solid histogram). 
The dashed curve and the dot-dashed curve represent 
the Maxwell distributions in three and two dimensions, 
which are defined by 

8 2 -/h 2 - 1)V2 

N »M=* 2Km + oLw ^-rt> (4) 
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Fig. 12. — Histogram of the kinetic energy of the downstream 
electrons within 430 < x < 442 at t = 233 (solid line). The 
horizontal axis is the kinetic energy normalized by the electron rest 
mass energy, 7 — 1. The dashed curve and the dot-dashed curve 
represent the Maxwell distribution for three and two dimensions, 
respectively. 
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Fig. 13. — Histogram of the kinetic energy of all the electrons 
in the simulation box (solid line). The three-dimensional Maxwell 
distribution is also plotted as a reference (dashed curve). There is 
no significant particle acceleration. 

and ^ 

N 2 n(7) = n d -^Le X p{-(~f-l)6], (5) 

respectively. Here, is the number density in the 
downstream region, Kj is the modified Bessel fun ction 
of the second kind (jAbramowitz fc Stegun I [l965t ) . and 
9 is determined so that 7^(7)^7 = (7) ~ \f5, 
which is expected if the upstream bulk kinetic energy 
is completely converted into the internal kinetic energy 
in the downstream rest frame. From Fig. [3^., we consider 
nd ~ 3.3 n e Q. It is observed that the electron distribution 
fits neither of the Maxwell distributions. 



Figure[13]shows the energy histogram of all electrons in 
the simulation box by a solid histogram. For reference, 
we plotted the three-dimensional Maxwell distribution 
by the dotted curve. Their normalizations are arbitrary, 
and there is no significant particle acceleration. 

4. CONCLUSIONS 

The simulation results clearly show that collisionlcss 
shocks can be driven in unmagnetized electron-positron 
plasmas. The structure of the collisionless shock prop- 
agates at an almost constant speed. The dissipation of 
the upstream bulk kinetic energy is mainly due to the 
deflection of particles by the strong magnetic field gen- 
erated around the shock transition region by the Weibel- 
type instability. It is remarkable that this shock has 
no corresponding linear wave and therefore it can be re- 
garded as a kind of "instability-driven" shock wave. We 
also observed the electric and magnetic field generation 
even in the upstream of the shock transition region due 
to the charge-separated current filaments together with 
the electrostatic beam instability. We found that a frac- 
tion of the downstream thermalized particles return up- 
stream through the shock transition region. These parti- 
cles would interact with the upstream incoming plasma 
and cause the generation of the charge-separated current 
filaments as well as the electrostatic beam instability up- 
stream of the transition region. We did not observe any 
efficient particle acceleration processes in our simulation, 
whereas we observed relatively rapid decay of the gen- 
crated magnetic field within the shock transition region. 
On the other hand, in the real world, there are many evi- 
dences to suggest that particles are accelerated at shocks, 
and there are several indications that large-scale, sub- 
equipartition magnetic fields may be generated around 
these shocks. The reasons why these processes were not 
observed in our simulation may be that (1) the back- 
ground magnetic field is essential for them, (2) the spa- 
tial and temporal scales of the simulation are too small 
to deal with them, or (3) these processes are inefficient 
in the electron-positron shocks. 
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